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Report  Developed  under  Small  Business  Innovation  Research  Contract.  A  Personal  Air  Ventilation  System  (PAVS) 
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Computational  fluid  dynamics  (CFD)  iterations  of  the  blower  were  used  to  minimize  losses  and  increase  efficiency 
through  adjustments  of  the  specific  blower  geometry.  CFD  outputs  included  the  total  and  static  pressure  rise,  as  well  as 
shaft  power  input.  The  shaft  power  requirement  was  matched  with  a  commercial  off-the-shelf  (COTS)  motor  such  that 
the  electrical  power  input  was  under  1 5  W.  The  motor  voltage  and  current  were  matched  with  standard  military  battery 
cells  of  various  chemistries,  yielding  run-times  that  exceeded  the  4-hr  run  time  requirement.  A  full  3D  system  model  was 
constructed,  and  the  final  system  volume  and  weight  (including  battery)  were  less  than  the  60  in3  and  2  lb  requirements. 
A  stereolithograph  mock-up  of  the  blower  was  built  and  delivered  to  demonstrate  functionality. 
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CENTRIFUGAL  BLOWER  FOR  THE 
PERSONAL  AIR  VENTILATION  SYSTEM  (PAVS)  -  PHASE  I 

1.  INTRODUCTION 


A  Personal  Air  Ventilation  System  (PAVS)  provides  cooling  to  an  individual  by  enhancing 
sweat  evaporation.  Even  though  the  ambient  air  may  be  warm  (i.e.  >  35°C),  provided  the  air  is 
relatively  dry,  there  is  significant  potential  to  remove  heat  as  illustrated  in  Figure  1. 


Blower  Heat  Removal  @95°F  (35°C) 


Exiting  Relative  Humidity 

Figure  1:  PAVS  heat  removal  potential  - 10  Cubic  Feet  per  Minute  (CFM),  95°F  ambient  air) 

Since  the  PAVS  system  is  extra  equipment  that  the  user  will  have  to  carry,  the  system  must  be 
compact  and  efficient  to  provide  a  true  benefit  to  the  user.  The  requirements  of  the  PAVS  in  this 
development  effort  are  summarized  as  follows: 

>  Air  Flow:  10  CFM  (0.28  m3/min)  through  (2)  C2A1  filters  (in  parallel) 

>  Backpressure:  5  inH20  (1243  Pa) 

>  Power  Consumption:  1 5  W  max 

>  Weight  &  Volume 

•  <2  lb,  <60  in3 

•  Includes  blower,  motor,  power  source 

•  Excludes  filters 

>  Depth  Dimension  <  2  in 

>  Run-time:  4  hr 

Given  the  15  W  maximum  power  consumption,  a  review  of  Commercial  Off-The-Shelf  (COTS) 
blowers  reveals  about  a  2X  factor  difference  between  available  performance  and  the  performance 
required  for  the  PAVS  (see  Figure  2).  Thus,  a  significant  efficiency  improvement  over  COTS 
performance  is  required. 
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2.  AXIAL-RADIAL  BLOWER 


A  blower  that  meets  the  requirements  of  this  development  effort  requires  a  high  aerodynamic 
efficiency.  The  feasibility  of  this  efficiency  can  be  examined  with  ns-ds  plots.  Ns  is  the  specific 
speed  of  the  turbo-machine,  and  ds  the  specific  diameter  (see  definitions  below).  The  plots  are 
estimated  maximum  achievable  efficiencies  based  on  historical  turbo-machinery  data.  It  should 
be  emphasized  that  they  do  NOT  necessarily  provide  a  precise  efficiency  for  a  given  turbo¬ 
machine,  but  they  do  give  a  “ball-park”  estimate  of  what  efficiencies  are  achievable. 


- —  specific  speed 

(H„gr 


,  p(g.,g)x 

VF 


specific  diameter 


co  is  angular  speed  (27i*rpm) 

V  is  volumetric  flow  rate 

Had  is  the  fan  pumping  head  (in  meters) 

g  is  the  gravitational  constant  (9.81  m/s  ) 


Figure  3  plots  the  maximum  efficiency  as  a  function  of  specific  speed  for  various  compressor 
types.  The  figure  indicates  that  for  the  specific  speed  of  the  proposed  blower  (indicated  with  red 
arrows),  a  “mixed  flow”  (mixed  axial/radial)  is  an  appropriate  choice.  Figure  4  is  a 
topographical  ns-ds  plot,  which  shows  that  the  maximum  achievable  aerodynamic  efficiency  for 
the  proposed  blower  should  be  ~  80%.  Thus  the  ns-ds  comparison  indicates  that  the  high 
aerodynamic  efficiency  required  for  the  proposed  blower  is  feasible. 
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Figure  3:  Maximum  compressor  efficiencies  as  a  function  of  specific  speed  (various  types) |1*. 


Figure  4:  nsds  diagram  for  single  stage  compressors  m. 
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3.  BLOWER  DESIGN 


The  design  of  the  blower  began  with  establishing  a  set  of  geometric  variables  used  to  construct 
the  blower.  These  variables  are  sub-divided  into  independent  and  dependent  categories.  A 
known  set  of  independent  variables  fully  constrains  the  geometric  design  of  the  blower.  In  order 
to  estimate  the  values  of  the  initial  set  of  independent  variables,  aerodynamic  similarity  of  turbo¬ 
machinery  was  used. 


Aerodynamic  Similarity 

The  proposed  blower  is  aerodynamically  similar  to  an  existing  centrifugal  compressor  pictured  in 
Figure  5.  The  performance  plot  of  this  compressor  demonstrates  a  high  aerodynamic  efficiency. 
Aerodynamic  similarity  of  turbo-machinery  is  measured  by  comparing  3  main  attributes: 

1)  Loading  Function  (\j /)  -  ratio  of  total  enthalpy  rise  (specific  work  transfer)  to  KE  of  blade 

2)  Flow  function  ((j))  -  flow  handling  capability  of  turbo-machine  (choked  flow) 

3)  Geometry  -  blade  height,  ID,  OD,  ratios 

The  loading  function  and  flow  function  are  defined  below.  The  3  aerodynamic  similarity 
attribute  values  are  highlighted 
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Table  1,  demonstrating  that  the  two  turbo-machines  are  similar. 

Mj, 

t/2 

Aho  =  change  in  stagnation  enthalpy 
U  =  blade  tip  speed 
m  =  mass  flow  rate 

y,  R  =  gas  constant  and  specific  heat  ratio  for  air 
Po  To  =  stagnation  temperature  and  pressure  at  inlet 
d  =  blade  tip  diameter  at  impeller  exit 


GAP  5% 
—  76300 
-•-82000 
-•-90250 
—•—99200 
-♦—103000 
— <  — 108300 


Inlet  guide 


vanes 


Impeller 


1.55- 

1.50- 

1.45  - 
1.40- 


0-  1.35  - 
CL 

1.30 


1.25- 


1.20- 


0.003  0  004  0.005  0.006  0.007  0.008  0.009  0.010  0  011  0.012  0.013 
m  *  sq  rt(RT  ^D'sqrtfr) 


diffuser 


Figure  5:  Photo  &  Performance  plot  of  existing  centrifugal  compressor 
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Table  1:  Comparison  of  existing  CC  and  proposed  blower 


Attribute 

Existing  CC 

Proposed  Blower 

Diameter  (cm) 

4.8 

4.4 

Blade  Height  (mm) 

2 

3.4 

Speed 

103,000 

25,000 

#  blades 

20 

10 

Pressure  Ratio 

1.46 

1.012 

Pressure  Rise  (Pa) 

46,000 

1243 

Mass  Flow  Rate  (g/s) 

6.8 

5.4 

Exit  Flowrate  (m3/min,  [CFM]) 

0.36,  ri31 

0.28,  F101 

Isentropic  Flow  Work  (W) 

278 

6.2 

Flow  function 

0.0073 

0.0069 

Isentropic  Efficiency 

73% 

55% 

Loading  function 

0.72 

0.63 

Motor  Efficiency 

92% 

80% 

Input  Power  (W) 

360 

14 

Geometric  Parameters 

Starting  with  values  established  from  aerodynamic  similarity,  a  set  of  independent  variables  was 
established  for  the  blower  as  shown  in  Table  2.  Also  listed  in  the  table  are  the  final  values  of  the 
blower  geometry  after  the  Computational  Fluid  Dynamics  (CFD)  analysis  and  design 
adjustments.  Many  of  the  variables  are  unchanged,  indicating  that  the  aerodynamic  similarity 
gave  a  reasonable  starting  point  for  the  blower  design. 

Table  2:  Independent  variables  (dimensions  in  mm,  angles  in  degrees  from  axial) 


Symbol 

Description 

Initial 

Estimate 

Final 

Value 

Impeller  inlet  hub  radius 

5 

8 

a; 

Impeller  inlet  hub  angle 

0 

0 

Impeller  inlet  blade  height  (radial) 

4 

4 

lb 

Impeller  length  (axial) 

10 

10 

ae 

Impeller  exit  hub  angle 

50 

50 

bc 

Impeller  exit  blade  height  (nonnal) 

2.5 

3 

Nb 

#  blades 

12 

12 

tb 

Blade  &  vane  thickness 

0.5 

0.5 

tr 

Inlet  radial  tip  gap 

0.08 

0.08 

t„ 

Exit  normal  tip  gap 

0.05 

0.05 

Nr 

Rotational  speed  (rpm) 

25,000 

35,000 

cw 

Radial  gap  b/n  rotor  &  stator  hubs 

1 

1 

fii.crat 

Diffuser  outside  radius 

25 

25 
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Symbol 

Description 

Initial 

Estimate 

Final 

Value 

by 

Vane  height 

2 

2 

Ids 

Straight  diffuser  length 

20 

20 

Ny 

#  vanes 

24 

24 

Ns 

#  Diffuser  vane  segments 

3 

3 

P2 

Relative  flow  angle  at  blade  exit 

20 

25 

Gt-sliu 

Blade  exit  to  vane  inlet  slip  angle  (relative) 

0 

6.3 

V 

Inlet  Volumetric  flowrate  (CFM) 

10 

10 

AP 

Blower  Pressure  Rise  (Pa) 

1250 

1250 

Ti 

Inlet  temperature 

35°C 

35°C 

P3 

Blower  Outlet  pressure 

1  atm 

1  atm 

A  few  of  the  important  dependent  variables  are  shown  in  Table  3  along  with  their  defining 
equations. 

Table  3:  Dependent  variables 


Symbol 

Description 

Equation 

1* 

Blade/vane  gap  length  (steamwise) 

3*be 

re 

Impeller  exit  hub  radius 

re  =  ri  +  y  [tan(a,.)  +  tan(ae)] 

ri,t 

Impeller  inlet  tip  radius 

bt 

r  =  r  H - - - 

cos(a;.) 

I*i,s 

Shroud  inlet  radius 

II 

+ 

re,t 

Impeller  exit  tip  radius 

re,t  =  re+becos(ae) 

t"e,s 

Shroud  exit  radius 

re,s  =re  +  (be  +  t„)cos(ae) 

lb,t 

Tip-to-tip  impeller  length 

kt  =lb-besm(ae) 

ls,t 

Shroud  length  (axial) 

ls,t=lb-(be+tn)Sin(ae) 

ldhi 

Axial  location  of  diffuser  hub  inlet 

1  -  I  +  °w 

1 dhi  lb  ^  .  ,  x 

tan  (ae) 

fitis 

Vane  inlet  shroud  radius 

rd,  =re,S+lgS'n(ae) 

Idti 

Vane  tip  inlet  shroud  axial  location 

ldti=ls,,+lgC°^ae) 

ldhi 

Diffuser  hub  inlet  radius 

rdln=re+CW 

t"dh2 

Vane  hub  inlet  radius 

rdh2=rdis~(re,s-re) 

ldli2 

Vane  hub  inlet  axial  location 

ldh2=Li+(be  +  t„)sHae) 

fdan 

Diffuser  inside  hub  radius  inlet 

r  —  r  —  h 

d,in  d  ,out  v 

1dm 

Axial  length  of  mixed  diffuser 

ldm=  dm  dhl  sm(ac) 
l-cos(aJ 
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Symbol 

Description 

Equation 

P-dm 

Radius  of  mixed  diffuser  arc 

D  _  ^ dm 

dm  ~  ■  ,  , 

sin  (ae) 

Pi 

Blower  Inlet  Pressure 

II 

1 

5 

4.  CFD  ANALYSIS 

Having  established  a  geometric  design  based  on  an  analytical  analysis  and  aerodynamic 
similarity,  the  next  step  is  to  estimate  the  blower  performance  using  a  Computational  Fluid 
Dynamics  (CFD)  program.  The  CFD  results  enabled  adjustment  of  geometric  constraints  to 
mitigate  losses  and  improve  efficiency. 

Surfaces  &  Meshing 

Since  the  blower  has  12  blades  on  the  impeller,  a  periodic  section  is  a  30°  angular  section. 
However,  the  30°  section  must  “twist”  to  follow  the  lines  of  symmetry  between  blades  and/or  cut 
through  the  middle  of  the  blades  and  vanes.  The  periodic  section  is  comprised  of  1  blade 
passage  and  2  vane  passages.  Figure  6  shows  a  wireframe  of  the  model  alongside  a  single 
periodic  section. 


Figure  6:  Wireframe  model  of  12  segments  (left),  and  a  single  periodic  section  (right) 

For  the  CFD  analysis,  only  the  flow  surfaces  are  important,  so  most  of  the  solid  material  can  be 
cut  away.  Additionally,  a  surface  must  be  added  for  the  analysis,  at  the  sliding  mesh  interface 
between  the  rotating  impeller  and  the  stationary  diffuser.  Various  surfaces  of  the  CFD  analysis 
are  shown  in  Figure  7. 
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Diffuser 


Figure  7:  CFD  defined  surfaces 


The  resolution  of  the  CFD  surface  mesh  needs  to  be  fine  enough  to  ensure  accuracy,  yet  coarse 
enough  to  run  a  simulation  in  a  reasonable  amount  of  time.  A  photo  of  the  initial  mesh  is  shown 
in  Figure  8. 


Simulation  Details 

The  analysis  involves  the  examination  of  the  unsteady  effects  due  to  flow  interaction  between  the 
stationary  components  and  the  rotating  blades.  In  this  study,  the  sliding  mesh  capability  of 
FLUENT  is  used  to  analyze  the  unsteady  flow  in  the  blower.  The  rotor-stator  interaction  is 
modeled  by  allowing  the  mesh  associated  with  the  rotor  blade  row  to  rotate  relative  to  the 
stationary  mesh  associated  with  the  stator  blade  row. 
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Various  boundary  conditions  are  placed  on  appropriate  surfaces  as  shown  in  Figure  9.  These 
include  specifying  the  mass  flow  rates,  outlet  pressure,  and  impeller  rpm  as  well  as  defining 
periodic  and  “walled”  surfaces.  Walled  surfaces  are  considered  to  be  insulated  with  a  “no-slip” 
condition. 


Inlet 


Some  information  specific  to  the  FLUENT  software  for  the  simulation  includes  the  following: 

a.  Density  Based  Solver 

b.  Implicit  Formulation 

c.  Unsteady 

d.  2nd-Order  Implicit  Unsteady  Formulation 

e.  Turbulence  model:  Renormalized  k-s 

f.  Near  wall  treatment:  wall  functions 

The  periodic  time  step  is  the  time  it  takes  for  consecutive  blades  of  the  impeller  to  pass  a  fixed 
node  on  the  diffuser.  This  time  step  is  a  function  of  the  rpm  and  #  blades  of  the  impeller: 

60 

Tperioi  ~  Nr  ■  Nb 

At  25,000  rpm  with  12  blades,  this  time  step  is  0.2  ms.  In  the  CFD  analysis,  the  initial  time  step 
is  an  integer  factor  smaller  than  this  number  (i.e.  1/4).  As  periodic  waveforms  begin  to  appear  in 
the  pressure  and  velocity  history  plots,  the  time  step  is  reduced  to  l/20th  this  number  for  final 
convergence  of  the  model.  Convergence  is  checked  by  comparing  2  consecutive  100  step 
calculations  for  similarity. 

The  CFD  analysis  generates  a  solution  by  solving  the  flow,  turbulence  and  energy  equations. 

The  air  is  assumed  to  be  an  ideal  gas,  and  the  SIMPLE  algorithm  (Semi-Implicit  Method  for 
Pressure-Linked  Equations)  is  used  for  pressure-velocity  coupling.  Second  order  up-wind  is 
used  for  density,  momentum,  turbulent  dissipation  rate,  turbulent  kinetic  energy  and  energy.  For 
pressure,  the  standard  discretization  scheme  is  used. 
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5.  CFD  RESULTS 


The  1st  CFD  run  was  completed  twice  with  2  meshes:  the  2nd  being  finer  than  the  1st.  Agreement 
between  the  two  meshes  suggested  that  the  mesh  resolution  was  sufficient  for  accurate  results. 
Results  of  the  CFD  analysis  are  analyzed  by  creating  the  following: 

1)  Tabulated  Velocity,  Static  Temperature  and  Pressure  (Ps,  Ts)  vs.  position  @  points  “0” 
(circular  entrance),  “1”  (impeller  entrance),  “2r”  (sliding  vane  mesh,  rotor  side),  “2s” 
(sliding  vane  mesh  diffuser  {stator}  side)  “3”  (diffuser  exit),  and  “4”  (circular  exit) 

2)  History  plots  of  velocity  and  temperature  vs.  time  step 

3)  Contour  and  streamline  plots  of  pressure  and  velocity 

Due  to  the  low  pressure  ratio  of  the  blower,  the  pressures  converge  faster  than  the  total  enthalpy 
used  to  calculate  the  required  shaft  power.  Thus,  when  it  became  apparent  in  the  1st  CFD  run  at 
25,000  rpm  that  the  pressure  rise  was  insufficient;  the  speed  was  increased  to  35,000  -  the  max 
speed  of  the  desired  motor.  This  was  the  2nd  CFD  run.  This  run  was  allowed  to  get  closer  to 
convergence  on  total  enthalpy  to  give  an  initial  shaft  power  estimate  for  the  blower. 

Each  successive  CFD  run  uses  the  same  conditions.  The  differences  between  runs  lie  in  the 
specific  geometry  such  as  blade/vane  angles  and  positions.  The  major  parameters  that  set  these 
values  are  summarized  in  Table  4.  The  table  also  contains  3  important  CFD  output  results:  gross 
total  pressure  rise,  net  static  pressure  rise,  and  shaft  power. 


Table  4:  CFD  Parameters  and  Results  for  Each  Run 


1st  Run 

2nd  Run 

3rd  Run 

4th  Run 

5th  Run 

6th  Run 

Target 

Impeller  rpm 

25,000 

35,000 

35,000 

35,000 

35,000 

35,000 

- 

r;  (mm) 

5.0 

5.0 

8.0 

8.0 

8.0 

8.0 

- 

P2  (deg) 

20 

20 

25 

25 

25 

25 

- 

be  (mm) 

2.5 

2.5 

3.0 

3.0 

3.0 

3.0 

- 

asiiD  (deg) 

0.0 

0.0 

0.0 

7.1 

6.0 

6.3 

- 

Gross  Total  AP  (Pa) 

495 

1199 

1884 

1896 

2011 

2013 

- 

Net  Static  AP  (Pa) 

170 

660 

240 

1108 

1243 

1253 

1243 

Shaft  Power  (W) 

n/a 

7.6 

10.8 

9.9 

10.2 

10.2 

11 

The  gross  total  pressure  rise  is  the  change  in  total  pressure  of  the  air  after  the  work  done  by  the 
impeller.  In  a  100%  efficient  blower  decelerating  the  flow  all  the  way  to  zero  velocity,  this  is  the 
maximum  possible  static  pressure  that  can  be  generated.  So  obviously,  this  number  should  be 
significantly  higher  than  the  1243  Pa  required  static  pressure  rise.  The  3rd  CFD  run  was  the  first 
run  where  the  total  pressure  rise  was  sufficiently  high,  but  the  losses  were  also  high,  and  the  net 
static  pressure  rise  was  only  240  Pa.  Thus,  in  the  4th  run,  asnp  (defined  as  the  relative  angle 
between  the  impeller  blade  exit  and  diffuser  vane  inlet)  was  adjusted  to  reduce  losses.  Runs  5 
and  6  further  adjusted  the  flow  angle  to  an  optimal  value. 

Detailed  tabulated  results  of  Run  6  are  shown  Table  5  at  various  plane  locations  of  the  blower. 
The  total  pressure  losses  (APt)  within  the  diffuser  are  a  measurement  of  flow  separation  and 
isentropic  inefficiency.  The  required  shaft  power  can  be  calculated  from  the  change  in  total 
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enthalpy  (Aht)  multiplied  by  the  mass  flow  rate  (5.3  g/s).  Because  only  the  impeller  does  work 
on  the  fluid,  the  only  change  in  total  enthalpy  should  be  from  1  to  2r.  The  non-zero  numbers  in 
the  table  reflect  the  accuracy  limits  of  the  analysis,  but  it  is  important  to  note  that  they  are  less 
than  2%  of  the  total  enthalpy  change. 


Table  5:  Tabulated  CFD  Results  (6th  Run) 


Plane 

ht 

Aht 

hs 

Ahs 

Pt 

APt 

Ps 

APs 

Total  temp 

ATt 

Static  temp 

ATs 

Velocity 

magnitude 

Axial 

Velocity 

.i/kg 

j/kg 

.i/kg 

.i/kg 

Pa 

Pa 

Pa 

Pa 

K 

K 

K 

K 

m/s 

m/s 

0 

circular  inlet  (0) 

9914.5 

9861.1 

100173.3 

100113 

308 

307.95 

11.18 

10.009 

1 

impeller  inlet  (1) 

9917.96 

3 

9640.5 

-221 

100168.4 

-5 

99853.9 

-259 

308 

0 

307.73 

-0.22 

22.495 

21.096 

2 

sliding  interface  rotor  (2r) 

11800.5 

1883 

11067.4 

1427 

102181.5 

2013 

101237.6 

1384 

309.9 

1.9 

309.15 

1.42 

39.346 

9.24 

3 

sliding  interface  stator  (2s' 

11800.3 

0 

11067 

0 

102179.5 

-2 

101238.5 

1 

309.9 

0 

309.15 

0 

39.65 

9.004 

4 

diffuser  vane  1  exit  (vl) 

11784.9 

-15 

11166.2 

99 

101759.3 

-420 

101053.2 

-185 

309.9 

0 

309.25 

0.1 

34.672 

18.363 

5 

diffuser  vane  3  exit  (3) 

11811.1 

26 

11620.3 

454 

101476.7 

-283 

101259.1 

206 

309.9 

0 

309.7 

0.45 

18.892 

18.284 

6 

circular  exit  (4) 

11825.1 

14 

11768.2 

148 

101427.6 

-49 

101366.1 

107 

309.9 

0 

309.83 

0.13 

10.54 

8.256 

Exit  -  Inlet 

1911 

1907 

1254 

1253 

1.9 

1.88 

Convergence  of  each  CFD  run  can  be  verified  by  examining  the  time  history  plots  of  velocity 
and  temperature  as  shown  in  Figure  10.  Note  that  the  temperature  only  begins  to  converge  after 
the  velocity  has  reached  near-convergence. 


Figure  10:  Convergence  Plot  of  Velocity  (left)  and  Temperature  (right) 
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Flow  angle  correction 


An  important  parameter  adjusted  in  subsequent  CFD  run  is  asup,  defined  as  the  relative  angle 
between  the  impeller  blade  exit  and  diffuser  vane  inlet.  The  adjustments  made  to  asup  are 
estimated  by  examining  the  velocity  vector  field  and  pressure  plot  outputs  of  the  analysis.  The 
adjustment  is  proportional  to  the  angle  measured  between  the  vane  and  the  velocity  vectors  when 
viewing  normal  to  the  flow  surface.  The  effect  of  increasing  asiip  is  to  “turn”  the  vane  into  the 
flow.  When  ashp  is  under-corrected,  a  high  pressure  stagnation  region  forms  “below”  the  vane, 
and  a  low-pressure  separation  region  occurs  “above”  the  vane.  When  asup  is  over-corrected,  the 
high  and  low  pressure  regions  flip-flop  to  opposite  sides  of  the  vane.  This  is  shown  pictorially  in 
Figure  11.  Figure  12  and  Figure  13  show  the  pressure  plots  of  Run  3  and  Run  4,  respectively. 
Note  that  Run  3  has  an  under-corrected  asup,  while  Run  4  is  slightly  over-corrected.  Figure  14 
shows  how  the  correction  angle  is  measured  in  Run  3  to  estimate  a  new  value  for  asiiP. 
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Contours  of  Static  Pressure  Ipascall  (Time=l  ,33l3e  02)  Sep  03,  2007 

FLUENT  6.3  13d,  dp,  pbns,  rke,  unsteady) 


Figure  12:  Pressure  Plot  under-corrected  asiip  (Run  #3) 


Contours  of  Static  Pressure  (pascal)  (Time  2.4 156e-02)  Sep  12.  2007 

FLUENT  6.3  (3d.  dp.  pbns.  rle.  unsteady) 

Figure  13:  Pressure  Plot,  over-corrected  aS|ip  (Run#4) 
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Velocity  Vectors  Colored  By  Velocity  Magnitude  (m/s)  (Time  =  1 .48 1  4eS@p)  03,  2007 

FLUENT  B.3  (3d,  dp,  pbns,  rke,  unsteady) 


7 . 0  6  e  + 
6.7 1  e  + 
6.36e  + 

6. 0  0  e  + 
5.65e  + 
5.30e  +  0 1 
4.95e  +  0 1 
4.59e  +  0 1 
4.24e  +  0 1 
3.89e  +  0 1 
3.54e  +  0 1 

3. 1  Be  +  0 1 
2.83e  +  0 1 
2.4  8e  +  0 1 

2.1  2e  +  0 1 
1 ,77e  +  0 1 
1 ,42e  +  0 1 
1 . 0  7e  +  0 1 

7. 1  4  e  +  0  0 
3.6 1  e  +  0  0 
8 . 3  8  e  -  0  2 


Vane 


Figure  14:  Flow  Separation  and  asiiP  correction  measurement  (Run  #3) 


Further  Blower  Improvements 

With  adjustment  and  refinement  of  the  blower  flow  angles,  the  feasibility  of  the  blower  design 
has  been  established.  While  there  is  still  room  for  tweaking  of  the  diffuser  design  (#  vane 
segments,  individual  vane  angles,  etc.),  the  next  logical  step  in  the  blower  development  is  to 
build  and  test  a  prototype  blower.  Test  results  will  then  be  compared  to  CFD  outputs  to  ensure 
agreement,  and  may  lead  to  adjustments  in  the  CFD  model.  Once  the  CFD  and  prototype  yield 
consistent  results,  then  final  parameter  tweaking  on  the  CFD  model  can  be  accomplished, 
resulting  in  a  2nd  improved  prototype  to  be  built  and  tested. 


6.  MOTOR  AND  BATTERIES 

The  Phase  I  design  shaft  power  requirement  is  10.2  W  at  35,000  rpm.  A  brushless  DC  motor 
from  Maxon  Motors  (#  201162)  has  been  identified  as  a  good  candidate  for  the  blower.  It  is  22 
mm  in  diameter  and  50  mm  long,  features  integrated  electronics,  and  weighs  only  85  g.  A 
technical  drawing  of  the  motor  is  shown  in  Figure  15.  This  motor  is  able  to  meet  the  power 
requirements  at  24  V  and  76%  efficiency  as  shown  in  Figure  16.  The  blower  has  been  designed 
such  that  the  impeller  will  mount  directly  onto  the  shaft  of  this  motor  with  no  additional  bearings 
required.  Thus,  the  electrical  power  required  is  10.2  /  0.76  =  13.5  W.  With  the  known  voltage 
power,  and  run  time,  a  suitable  battery  can  be  identified. 
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OO  OperPt. 


Torque  [mN*m] 

-  24  V 

OO  OperPt. 

Figure  16:  Motor  Performance  at  24  V  with  operating  point  indicated 


Table  6  lists  a  wide  range  of  suitable  battery  cell  combinations  based  on  4  chemistries  available 
from  Saft®:  Li-MnCL,  Li-SiCL,  Li-SOCh,  and  rechargeable  Li-ion.  This  battery  analysis  uses 
just  the  cells  available  from  one  manufacturer,  but  is  representative  of  what  is  available  in  the 
market.  The  table  lists  only  those  combinations  that  can  meet  the  24  voltage  requirement  and 
run  for  at  least  4  hours.  The  “best”  cell  combination  for  each  battery  chemistry  is  highlighted  in 
Table  6.  The  Li-SOCL  has  by  far  the  highest  energy  density,  and  would  last  over  10  hours,  but 
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is  also  the  most  toxic.  The  Li-ion  combination  lasts  over  5  hours  and  has  the  advantage  of  being 
rechargeable,  but  with  the  penalty  of  increased  weight.  Moving  forward  in  the  current  size  & 
weight  analysis,  a  representative  battery  will  be  assumed  to  weigh  400  g  and  have  the  form 
factor  of  9  “C”  sized  cells. 

Table  6:  Suitable  Saft®  battery  cell  combinations 

Nominal  Max  Nominal  Nominal  Cell  Energy  Nominal  Total 


Chemistry 

Cell  Type 

Cell  Size 

Nominal 
Voltage  (V) 

Drain 

(mA) 

Current 

(A) 

Capacity 

(Ah) 

Energy 

(W-hr) 

Weight 

(9) 

Density 

(W-hr/kg) 

Battery 

Configuration 

Battery 

Voltage 

Battery 

Current 

Drain 

/cell 

Run  Time  Weight 
(hr)  (gram) 

Li-Mn02 

LM  17130 

1/3  A 

2.7 

4.5 

0.3 

0.5 

1.35 

8 

169 

9 

S 

5 

P 

24.3 

0.56 

0.11 

4.5 

360 

Li-Mn02 

LM  22150 

1/3  sub-C 

2.7 

40 

0.4 

0.9 

2.43 

15 

162 

9 

S 

3 

P 

24.3 

0.56 

0.19 

4.9 

405 

Li-Mn02 

LM  26500 

C 

2.7 

200 

1.5 

4.5 

12.15 

60 

203 

9 

s 

1 

P 

24.3 

0.56 

0.56 

8.1 

540 

l_i-S02 

G  04/3 

1/2  AA 

2.8 

50 

0.25 

0.45 

1.26 

8 

157.5 

9 

s 

5 

P 

25.2 

0.54 

0.11 

4.2 

360 

Li-S02 

G  06/2 

AA 

2.8 

80 

0.5 

0.95 

2.66 

15 

177 

9 

s 

3 

P 

25.2 

0.54 

0.18 

5.3 

405 

l_i-S02 

G  32/3 

2/3  A 

2.8 

80 

0.75 

0.8 

2.24 

12 

187 

9 

s 

3 

P 

25.2 

0.54 

0.18 

4.5 

324 

Li-S02 

G  36/2 

"long"  A 

2.8 

80 

1.5 

1.7 

4.76 

18 

264 

9 

s 

2 

P 

25.2 

0.54 

0.27 

6.3 

324 

l_i-S02 

LO  34  SX 

1/3  C 

2.8 

80 

1.0 

0.86 

2.408 

18 

134 

9 

s 

3 

P 

25.2 

0.54 

0.18 

4.8 

486 

Li-S02 

G  52/3 

C 

2.8 

1000 

2.5 

3.2 

8.96 

47 

191 

9 

s 

1 

P 

25.2 

0.54 

0.54 

6.0 

423 

l_i-S02 

LO  29  SHX 

C 

2.8 

250 

2.5 

3.75 

10.5 

40 

263 

9 

s 

1 

P 

25.2 

0.54 

0.54 

7.0 

360 

Li-S02 

G  54/3 

5/4  C 

2.8 

200 

2.5 

5 

14 

58 

241 

9 

s 

1 

P 

25.2 

0.54 

0.54 

9.3 

522 

l_i-S02 

LO  43  SHX 

5/4  C 

2.8 

200 

2.5 

5 

14 

53 

264 

9 

s 

1 

P 

25.2 

0.54 

0.54 

9.3 

477 

l_i-S02 

LO  40  SX 

2/3  "Thin"  D 

2.8 

120 

2 

3.5 

9.8 

40 

245 

9 

s 

1 

P 

25.2 

0.54 

0.54 

6.5 

360 

Li-SOCI2 

LS  26180 

1/3  C 

3.6 

10 

0.4 

1.2 

4.32 

24 

180 

7 

s 

2 

P 

25.2 

0.54 

0.27 

4.5 

336 

Li-SOCI2 

LSH  14 

C 

3.6 

15 

1.3 

5.5 

19.8 

51 

388 

7 

s 

1 

P 

25.2 

0.54 

0.54 

10.3 

357 

Li-ion  (rechargeable) 

MP  144350 

14.5x43x50mm 

3.75 

500 

2.6 

2.6 

9.75 

68 

143 

7 

s 

1 

P 

26.25 

0.51 

0.51 

5.1 

476 

7.  BLOWER  SYSTEM  MODEL 

A  3D  model  of  the  complete  blower  design  has  been  made,  including  the  battery  cells.  Figure  17 
shows  pictures  and  cross-sections  of  the  model.  A  weight  and  volume  estimate  of  the  blower  is 
shown  in  Table  7.  The  volume  is  the  envelope  of  space  occupied  by  the  blower,  such  that  a  2” 
diameter  is  assumed  to  occupy  a  2”  x  2”  square.  Since  the  motor  is  contained  inside  the  blower, 
it  does  not  contribute  to  the  volume. 
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Figure  17:  Final  Phase  I  Design  blower  model 


Table  7:  Blower  Weight  and  Volume  Estimate 


Material 

Weight 

Envelope 

Volume 

g 

lb 

(in3) 

Blower 

Aluminum 

268 

0.59 

19 

Motor 

Actual 

85 

0.19 

- 

Filter  Adapter 

Nylon 

42 

0.09 

9 

Battery 

Cells 

400 

0.88 

18 

TOTAL 

- 

795 

1.75 

46 

18 


8.  CONCLUSIONS 


This  Phase  1  effort  has  successfully  established  the  feasibility  of  the  PAVS  blower.  Table  8 
shows  that  the  Phase  I  design  meets  or  exceeds  all  of  the  weight,  volume,  and  performance 
requirements.  A  stereo-lithograph  model  of  the  blower  was  built  and  delivered,  demonstrating 
basic  functionality.  At  this  point,  a  prototype  blower  has  been  designed  and  is  ready  to  be  built 
and  tested  as  part  of  a  Phase  II  effort. 


Table  8:  Phase  I  Design  vs.  Requirements 


Parameter 

Requirement 

Phase  I  Design 

Air  Flow  (CFM) 

10 

10 

Pressure  Rise  (Pa) 

1243 

1253 

Electrical  Power  (W) 

15 

13.5 

Run  Time  (hr) 

4 

5  to  10 

(depends  on  battery  chemistry) 

Weight  (lb) 

2 

1.8 

Volume  (in3) 

60 

46 

Depth  (in) 

2 

2 
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